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Abstract 

In this paper we eompute and diseuss the differential eross seetion of the Bhabha seattering 
in the framework of the z = 2 Lifshitz QED. We start by eonstrueting the elassieal solutions 
for the fermionie fields, in partieular the eompleteness relations, and also deriving the theory’s 
propagators. Afterwards, we eompute the photon exehange and pair annihilation eontributions 
for the Bhabha’s proeess, and upon the results we establish the magnitude of the theory’s free 
parameter by looking for small deviations of the QED tree results. 
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1 Introduction 


Within the field theory eontext the Lifshitz models, based on an anisotropy between spaee and 
time Q, i.e. the sealing x‘ —)■ Ax‘ whereas t — KH, always had a variety of applieations permeating 
the literature of eondensate matter physies, and nowadays has been reviver and gained much attention 
within the context of the search for a perturbatively consistent theory for gravity Q. The remarkable 
advantage of this approach is mainly motivated by the possibility of defining new renormalizable 
interactions, but they are constructed in such a ways to avoids the appearance of ghosts (negative 
energy modes) whose presence is an imprint of theories with higher time derivatives. 

Although these Lifshitz field theories have a better ultraviolet behavior at the expenses of breaking 
Lorentz invariance, many interesting and rich features regarding the theory’s renormalizability [[^ and 
also in discussions on Lorentz symmetry violation 0 have been observed. In exploring implications 
of the time-space asymmetry of a given field theory’s renormalizability, many investigations were 
undertaken and reported recently, for instance, studies considering Lifshitz gauge field theories [j5]-|7|, 
in particular a detailed study of its general properties can be found in the Ref. 0Q; moreover, scalar 
field theories [[S -10|, four-fermion theory [ 1 L 121, and model 1131 were also considered. In 

completion, we may also cite the study of formal aspects of Lifshitz field theories according to the 
BPHZ scheme [ |T4l . 

Recently, an interesting study of a phenomenological viable z = 2 Lifshitz modified quantum 
electrodynamics (QED) has been reported p3| . By the fact that the theory is super-renormalizable, 
non-trivial aspects of its quantum content were discussed. Nevertheless, the possible effects on scat¬ 
tering processes in the framework of QED by these new super-renormalizable couplings have not 
been investigated. Some attention has been paid in studying scattering process, in particular to the 
Bhabha scattering, in the Eorentz violating field theories, due to the emergence of new vertices in the 
theory [ 16-181, and also at the generalized electrodynamics [ 191. In these studies, experimental data 
was used in order to obtain a bound value for the theories free parameters. 

Remarkably, Bhabha scattering is one of the most fundamental reactions in QED processes [20| 
and has played an important role in experiments in colliders [[2T]-(^, since it is particularly used 


to determine the luminosity of the collisions [24|. It is important to emphasize that nowadays there 


are ongoing theoretical analysis present in the literature [25 -271, in the search of an improvement in 
the calculations of higher-order contributions for the Bhabha’s cross section. Therefore, motived by 
theoretical and experimental grounds, we believe to be interesting and rather natural to consider and 
analyze the behavior of the differential cross section for Bhabha scattering in the framework of the 
z = 2 Eifshitz QED. 

In this paper, we discuss a scattering process in the framework of z = 2 Eifshitz modified quantum 
electrodynamics. In Sec. we start by making a brief review of the z = 2 QED and then we construct 
in details the classical solutions for the fermionic fields, as well as the completeness relations for 
the spinors, and also determine the theory’s fermionic and gauge field propagators quantities that 
are all needed in the calculation of the cross-section for the Bhabha scattering. Next, in Sec. 
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we discuss and present a detailed calculation of the S-matrix elements, the transition amplitude for 
the theory. The calculation is performed by taking into account the high energy limit, i.e. we take 
the fermionic mass to be zero; moreover, we express the kinematic variables in the center-of-mass 
frame. In conclusion, we establish the magnitude of the theory’s free parameter by looking for small 
deviations of the QED tree results [241. In Sec. |^we summarize the results, and present our final 
remarks and prospects. 


2 General discussion 


In this section we will introduce our basic notation and describe the model and its classical solu¬ 


tions. Let us consider here the following Lagrangian density for the z — l Lifshitz QED [ 151, 


1 


1 


C, = -Fo/Eo,' - -Fij {M^ - A) Fij + t/r {ijoDo - iMy^Dk - D^Dk - nF) Xf/, 


( 2 . 1 ) 


with the metric signature (-h, the covariant derivative 0^ = 9^ + igA^, and the field stress 

tensor F^y = d^Av-dyA^. This theory is known to be super-renormalizable, but still contains power 


counting diverging graphs, which were shown to be actually UV finite [15|. Erom the Lagrangian 


density (|2.1|) we may read that the length dimensions are, in 3 -1- 1 dimensions. 


[Ao] = [w]=L 2 , [Ai]=L -2, 
[g]=L-2, [m] = [M]=L-K 


The theory is invariant under the t/ (1) gauge symmetry 


1 


Af,^A^ + -d^a. 
§ 


( 2 . 2 ) 


(2.3) 


We are interested here in evaluate the Bhabha’s differential cross-section, and for that matter we need 
to construct the solution for the fermionic fields [ |2^ . Hence, in order to construct these solutions for 
the Dirac’s fields we shall consider the following modified free field equation, 

(iyo^o - iMjkdk - dkdk -m^)y/ = 0, (2.4) 

the coefficients of the differential equation are constants, thus V^(x) = (p ); will be a solution. 

Hence, we obtain 

{YoPo - MYkPk + P^-m^)x{p)= 0, (2.5) 

where = piPi. Eor the energy eigenvalues, we have 


Pq = ±Ep = ±\J M^p'^ F (p2 _ m2) , 


( 2 . 6 ) 


for each value of pQ, the solution ( |2.5| ) has a two-dimensional solution space. Eor the explicit calcu¬ 
lation, we choose 

Us {p) 


x{p) = 


Vs {p) 


(2.7) 








and using the following representation for the y-matrices, 





0 a A 
-do)' 


where 1 is a two-dimensional identity matrix, and d are the set of Pauli matrices. These all lead to, 


{Po +P^- m^) Us {p) = M (d.p) (p) , (2.8) 

{po -P^ + m^) Vi (p) = M (d.p) Us (p) . (2.9) 


Nonetheless, it is not difficult to show that u (p) corresponds to those solutions with positive energy, 
Po = -\-Ep, while V {p) are solutions with negative energy, po = —Ep. The next step involves in 
defining the energy projection operators, in such a way 


A+ (p) Us (p) = Us (p) , (2.10) 

(p) Vi (p) = Vi (p), (2.11) 


and 

(p) Us (p) = A+ (p) Vs (p) = 0. 

We can then show that the operators 

._TYoPo^MYkPk + P^-m^ 

2(p2_m2) 

satisfy the above relations, as well as 

A+(p)+A^(p) = 1 , A+(p)A^(p) =0, 

[A±(p)]- = A±(p). 


Thus, we have the completeness relations satisfied by 


Kp (^) = E (P’ (P’ = 

5=1 

2 

^ap (P) = -E^a (pa)v/3(p, 


21 


S = 


-YoPo + MYkPk + p^-m 

2 (p2 — ^2) 

7opo - MYkPk + p^-m^ 


2(p2 — m2) 


ap 

J ap 


Finally, we have the free solutions are written as 


= E 

r ' 

=E 


d^p 

(2;r)2 


p^ — trd^ ^ 


d^p f p^ — trd^ ^ 


(2;r)2 


br{p)Ur{p)e-^P^ + dl{p)Vr{p)e^P^ , 
bl (p) Ur (p) -h dr (p) Vr (p) , 


( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 
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where, the operators algebra read, 

\br{p)M (^)l = \dr{p),4 


= 5rsS (p-q). 


( 2 . 20 ) 


From sueh construction one can easily show that the equal-time anti-commutation relations is satisfied 


^ J xo=yo 


Besides, we can also determine the fermion propagator 


Sipo^p) =i 


.YoPo-MYkPk + P^ - trr 


Pq — — (p2 — frp-) 


We can proceed similarly for the gauge field, and find for the energy eigenvalues, 

kQ = itOyt = i y/ 

and the solutions 

d^k ( 1 


(^) = E 


(2^)- 


2Wi 


1 

2 r 


ai{k)£l{k)e '^^ + al{k)£l{—k)e 


ikx 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


(2.24) 


However, in order to define the photon propagator, we must impose a gauge condition, which reads a 
Lorentz-like condition 

Q. [A] = <9oAo - (-A-fM^) dkAk = 0, 


this leads to a non-local gauge-fixing term in the Lagrangian density [29|, = 1), 

1 


Cgf = - (^oAo - (-A-fM ) dkAk) 


2(-A + M2) 

but to a well-behaved expression for the propagator: 

{ko^k) = Dqq {ko, k) = 


{dQAQ-{-A-YM^)djAj), 


M^ + k^ 


k^-M^k^-k^' 


iT>ij (ko, k) = Dij (ko, k) = 


% 


kl-M^k^-k^' 


(2.25) 

(2.26) 
(2.27) 


with no off-diagonal components. Nevertheless, it was shown in the Ref. [ 181 that off-diagonal com¬ 
ponents are generated after quantum corrections. We can read from the Lagrangian density ( |2.1[ ) that 
the z = 2 QED presents three vertices: two three-point vertices. 


V/AoV/ ^ igro, wAkW -ig {jkM -f Ipf''^ + pf ^ 

which can be conveniently rewritten as in a compact form, 

(WAqY, WAk¥) Aa id) = ig (7o, -YkM - qk) , 


(2.28) 


(2.29) 


where we have introduced, by means of convenience in notation, the index a = 0,..., 3, which should 
not be confused with spacetime index, and qk = -l- and one four-point vertex, xftxi/AiAj — 

—ligdij. But as we are here interested in the Bhabha scattering, we shall need to consider only ( |2.29| ). 
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(a) 


(b) 


Figure 1: Bhabha scattering: (a) photon exchange and (b) pair annihilation process. 


3 Bhabha scattering 

Collision experiments in high energy physics provide has been served as a ground framework for 
testing many theoretical proposals in QFT. Moreover, Bhabha scattering, —)■ is one of the 


most fundamental reactions in QED processes [20|, and its experimental data has been extensively 


used in order to obtain bounds in the context of Lorentz violating models [ 16-18| as well as in the 
generalized electrodynamics eg- 

Let us consider the following description as for our system depicted at the Figj^ the spin and 
momenta for the ingoing {p\) and outgoing {p 2 ) electrons are {ri,pi) and (r/,P 2 ), respectively, 
while the spin and momenta for the ingoing (^ 1 ) and outgoing (^ 2 ) positrons are {si,qi) and 
{sf,q 2 ), respectively. Nonetheless, it should be remarked that, because the photon propagator does 
not have off-diagonal components, we shall not have a mixed contribution between the temporal Aq 
and spatial vertices in our calculation of the differential cross-section for the Bhabha scattering. 
From that we may now evaluate the total contribution for the transition amplitude, to lowest order. 


that corresponds to sum of the scattering diagram (s-channel) 

iM' = -Ng^Vsi (<?l) Aa [pi - qi)Uri {p\)Dab {pi +qi)Urf (Pl) {P2 “ qi) Vsf (qi) , 

and of annihilation diagram (t-channel) 

iM^ = Ng^Vsf {qi)K (<?1 + qi) Vsf (qi) Dab (Pl - Pl) Mr/ (Pl) Afe (pi -h P2) Ur^ (pi) , 
where we have defined 


N = 


2 2 \ 7 

p\ — m ^ ^ 

^Pi 


2 2 \ 7 

P2 ~ Ml ' ^ 
^P2 


2 2 \ 7 

q\ — m^^ ^ 
p 


2 2X7 

q2 — m^^ ^ 

p 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


In order to evaluate the cross section for the Bhabha scattering, we shall now compute |iAdf by 
taking an average over the spin of the incoming particles and sum over the outgoing particles 

1 




(3.5) 


ri,Si rf,sf 
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In this way, it follows 
1 


11 = 7 1 l-^f+T I |^T+T I {MrM!+j I {My 


(3.6) 




spin 


spin 


spin 


spin 


Although the calculation of each one of these terms is straightforward, the details are lengthy, but 
can be accomplished by using the completeness relations ( |2.17[ ) and by evaluating the traces of Dirac 
matrices products. Furthermore, as our main goal is to consider the behavior of the scattering process 
in the high energy limit, we shall take = 0 in our calculation. Nevertheless, in order to complete 
the calculation is interesting to express the kinematic variables in the center-of-mass frame, in which 


Pi = {E,v), ?1 = (£,-P), P2 = (£,q), ^2 = (£,-q) 


(3.7) 


Besides, due energy conservation: Ei = Ej |p| = |q|, and also we define 6 as the center-of-mass 
scattering angle, i.e. (p.q) = p^cosO. Hence, the resulting expression for all these contributions read, 
for the s-channel contribution. 


11-^^ 


1 

4 

2 g 






spin 

4 


{AcpyM^ 

1 


/,(£,M,0), 


it{E,M,e) 


spin 


83^^ sin^ (f) [1 +2 (1 - A) sin^ (|)]^ (1 - A)^ 


i E (m-ym ;+i E (mYm- = 


1 


4(£,M,0) 


■spin spin 8M4sin2|(l+2(l-A)sin2|) (1-A)(4 (p)3’ 

where we have defined the quantities Is, It and 1st, in the following way 


(3.8) 

(3.9) 
(3.10) 


Is {E,M, 0) =23 + 104(p + 76(p2 - 23A - 58(pA + 4(pcos 0 (1 - (p) 

+ cos 20 (23 + 84(P -h 32(p^ - 23A - 38(pA), (3.11) 


while 


It {E,M, 0) =2 (63 [1 - A] + (p (296 - 170A+ (p [299 + 32(p - 72A])) 

+ cos 0 (135 [1 - A] -h (P (621 - 351A + (P [536 - 96A] )) 

+ cos 20 (61 [1 - A] + 2(p (154 - 93A + 8(p [23 -f 4(p - 7A])) 

+ cos30(5(l-A) + (p(19 + 8(p-9A))+cos40((l-A)-f2(p(2 + (p-A)). (3.12) 


and 


1st {E,M, 0) =142 (1 - A) + (p [747 - 463A -h 2(p (385 + 72(p - 16A)] 

+ 2cos 0 [154 (1 - A) + (p (1241 - 933A + 2(p [1063 + 72(p - 252A])] 

+ cos 20 [182 (1 - A) + (p (783 - 419A + 2(p [349 -f 72(p - lOOA])] 

+ 8 cos 30 [11 (1 - A) + (p (52 - 30A-f 47(p - 9(pA)] 

+ 4cos40 [3 (1 - A) + 2(p (7 - 4A + (p [6 - A])]. (3.13) 
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with (p = ^ and A = +4-(p. Nevertheless, we have that the differential eross seetion (in natural 

units) is given by 

^= 52 iIEl‘-^l , (3.14) 

then, with the previous results one obtain. 


<■ '■/>^/ 


da 


a 


— =^h{E,M,e) + 




4(£,M,0) 


4(j!) 


+ 


;i-A) <?>sin4(f) [1+2(1-A)sin2(f)]^ 




ist{E,M,e) 


32 (1 - A) (p2 sin2 I (1 + 2 (1 - A) sin^ f) ’ 


(3.15) 


pZ 

where oc = |^. To eonelude this seetion, based on experimental grounds, we may determine limit 
bounds for the parameter M in the ease previously evaluated. For Bhabha seattering, the experimental 
data on preeision tests in QED are readily available in Ref. [24|, and small deviations on the QED 
tree results may be expressed in the form 



-1 

QED 


3Et 


hi 


(3.16) 


■± 


where Ecm = 2E = 29 GeV and A+ = 200 GeV with 95% eonfidenee level; aetually, A and is a pa¬ 
rameter representing possible experimental departures from the theoretieal predietions. Considering 
the leading eorreetions from (|3.15|) with |eos 0| < 0.55 [|24l, we ean show that the magnitude of these 
eorreetions are of order (p = ^, and therefore when we eompare them in d we obtain the bound limit 
4/2 > lO^^eV, for A = 200 GeV. 


4 Concluding remarks 

In this paper we have diseussed the Bhabha seattering in the framework of Eifshitz field theory, 
the z = 2 quantum eleetrodynamies. This theory has shown to be rather interesting both in elassieal 
and quantum eontext, and also phenomenologieal viable Eifshitz extension of QED. Although it is a 
super-renormalizable field theory, it provided rieh outeomes in the radiative eorreetions. 

Reeently, due to the aeeuraey of the experimental data of Bhabha seattering, interesting stud¬ 
ies have used this data in a stringent way in setting bound values on theory’s free parameters have 
appeared in the literature, sueh as in the ease of Eorentz violating theories and generalized eleetro¬ 
dynamies. Nonetheless, we have started by diseussing the elassieal theory in order to eonstruet the 
free solution for the fermionie fields, sinee they now obey a modified Dirae equation and possess a 
different dispersion relation. These results were neeessary in order to eompute eorreetly the S-matrix 
elements (transition amplitude). In partieular, we were interested in the high-energy behavior of the 
eross seetion, henee we took m2 = 0 in our ealeulation. Afterwards, we have eompared the leading 
eorreetions from the z = 2 QED with small deviations from the experimental data of the Bhabha seat¬ 
tering 95% eonfidenee level, finding that for values A/2 > lO^^eV the theory is eompatible with the 
experimental data. 
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Due to the richness of outcomes that z = 2 QED has provided, we believe that may be there are 
many other issues that should be discussed. For instance, we may cite many investigations on Lifshitz- 
like theories concerning in the study of the effective potential considering an arbitrary value of the 
critical exponent z and for different kind of couplings and fields [3^36|. In particular, an extension 


of our studies for z = 2 QED for the finite-temperature case [37| certainly would be very important, 
since the high-temperature corrections provide a description of a large class of interesting phenomena 
and these contributions may be put in a closed form expression [ [38||^ . These issues and others will 
be further elaborated, investigated and reported elsewhere. 
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